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\jQ [ Abstract 
T— I ■ 

' The background field formalism is used to implement nonpertur- 

bative QCD contributions into diagrammatic technic at T > 0. The 
leading terms both in the confining and nonconfining phase are iden- 
Q^ ■ tified at large Nc and the transition temperature is calculated in this 

. limit, which appears to be in good agreement with the Monte Carlo 

Qh! calculations. 

1 Introduction 

^ . It is believed that the perturbative QCD is applicable in the deconfined phase 

^ I at large enough temperatures T, where the effective coupling constant g(T) 

is small [1], while at small T (in the confined phase) the nonperturbative 
effects instead are most important. However even at large T the physics is 
not that simple: some effects, like screening (electric gluon mass), need a 
resummation of the perturbative series [2], while the effects connected with 
the magnetic gluon mass demonstrate the infrared divergence of the series 
[3]. 

There is a substantial amount of data from lattice calculations in the 
deconfined phase, which can be explained only by nonperturbative effects 
(see discussion below in section 6). Therefore one can assume that non- 
perturbative effects are of vital importance both for small and large T and 
should be taken into account before any dynamical perturbative scheme is 
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applied. We suggest to use for this purpose the background field formalism 
[4], in connection with the cluster expansion for background fields [5] and 
the Feynman-Schwinger representation (FSR) [6] for the Green's function. 

This method has been used successfully for QCD at T = and it contains 
confinement as an essential ingredient, appearing through the area law of 
Wilson loops in FSR. 

It is a purpose of this paper to apply the method to the case of T > 
by proper modifications of FSR and the background field method (BFM). 
For this purpose we write in Section 2 the familiar BFM expressions for 
the free energy with the only important new element: we do not consider 
the background field as classical (as it is usually done) but rather make an 
average over an ensemble of background fields, so that the latter enter final 
expressions for the free energy in the form of gauge-invariant correlators. 

As a consequence the simple loop diagram- the leading term 0{g^) of the 
usual perturbation theory — becomes a set of diagrams representing multiple 
loops interacting via nonperturbative correlators. This is done in Section 3 
for gluons and in Section 4 for quarks. 

The important new element there is a new form of FSR for a gluon or 
a quark Green's function appropriate for T > 0: it contains a path inte- 
gral winding on a torus which compactifies the fourth Euclidean direction. 
This enables us to compute gauge-invariant Green's function for hadrons and 
response functions for T > 0. 

Higher order amplitudes are discussed in Section 5. Here the behaviour 
of the QCD coupling constant g in strong background field at all distances 
is exploited [7]. 

In Section 6 we identify the leading contributions to the free energy in the 
limit when Nc ^ oo and g is small. It is gratifying that the phase transition 
of the first order occurs already in this oversimplifying case and the resulting 
transition temperature Tc does not depend on Nc when Nc is large. 

Clearly the mechanism of the phase transition is the evaporation of a 
part of the gluonic condensate into a gas of almost free gluons [8]. The 
predicted value of Tc agrees well with Monte Carlo calculations both in the 
case of gluodynamics [9] and in the case of QCD [10] for different number of 
flavours. 

In conclusion possible prospectives of the proposed method are outlined. 
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2 Basic equations 

We start with standard formulas of the background field formalism [4] gener- 
alized to the case of nonzero temperature. We assume that the gluonic field 
An can be split into the background field 5^ and the quantum field 

Ai,^Bi, + a^, (1) 

both satisfying periodic boundary conditions 

B,^{zi, Zi) = Bi^{z4 + nP, Zi), a^izi, zi) = 0^(2:4 + Zi), (2) 

where n is an integer and /3 = 1/T. 
The partition function can be written as 

z(i^,r,// = o) =< z{B) >B, 

Z{B) = N j D(t)exp{- dr j d^xL{x,T)) (3) 

where denotes all set of fields a^, N is a. normalization constant, 

and the sign <>b means some averaging over (nonperturbative) background 
fields Bfj^, the exact form of this averaging is not needed for our purposes. 
Furthermore, we have 

8 
i=l 

where 

L3 = e''{D\B)Ue'; L, = -igQ''{Dn, an)abQ' (4) 
L5 = -a{D^,{B)af,f\ Lq = Lint{a^, a'^) 
Lr = -a,Dn{B)Fn,{By, Lg = *+(m + D{B + a))^ 

Here ©, © are ghost fields, a- gauge-fixing constant, Lq contains 3-gluon- 
and 4-gluon vertices, and we keep the most general background field B^, not 
satisfying classical equations, hence the the appearence of L7. 
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The inverse gluon propagator in the background gauge is 



(5) 



where 



{Dx)ca = dx5ca - igT^Bl = dx5,a - gfbcaBl 



(6) 



We cosider first the case of pure gluodinamics, = 0. 
Integration over ghost and gluon degrees of freedom in (3) yields 



Z{B)^N'idetW{B));,f[det{-D,{B)D,{B + a))l^^ x 

oo cl X 1 



(7) 



One can consider strong background fields, so that gB^j, is large (as com- 
pared to A.QQ^), while ccg = |^ in that strong background is small at all 
distances [7]. 

In this case Eq. (7) is a perturbative sum in powers of gf", arising from 
expansion in (gia^)". 

In what follows we shall discuss the Feynman graphs for the free energy 
F(T), connected to Z{B) via 



As will be seen, the lowest order graphs already contain a nontrivial 
dynamical mechanism for the deconfinement transition, and those will be 
considered in the next section. 

3 The lowest order gluon contribution 

To the lowest order in gu/j, ( and keeping all dependence on gBf^ explicitly) 
one has 



F{T) 



Tin < Z{B) >B 



(8) 



Zo = e-^°(^)/^ = TV' < exp(-Fo{B)/T) >b 
where using (7) Fq{B) can be written as 



(9) 



^Fo(B) = -IndetW - lndet{-D\B)) 




(10) 
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In (10) Sp implies summing over all variables (Lorentz and color indices 
and coordinates), ({t) = Hm-^ \s=o is a regularizing factor, one can use 
also the Pauli-Villars form for ((t). 

Graphically the first term on the r.h.s. of (10) is shown in Fig.l and is a 
gluon loop in the background field, the second term is a ghost loop, shown 
in Fig. 2. 

Let us turn now to the averaging procedure in (9). With the notation 
(p — —Fq[B)/T, one can exploit in (9) the cluster expansion [5] 

oo 2 

< exp(£i > B= exp —7 = 

exp{< ^ >B +^[< ^''>B-<V >V\ + 0((^=^)}. (11) 

To get a closer look at < (^9 >b we first should discuss tcrmal propagators of 
gluon and ghost in the background field. We start with the thermal ghost 
propagator and write the FSR for it [6] 

{-DXl X I dte*^^(^) \y>^l dtiDz)le-^^x,y) (12) 

Here $ is the parallel transporter in the adjoint representation along the 
trajectory of the ghost: 

l>(x,y) = Pexp{ig J B^{z)dz^), (13) 

also 

4 70 ^' ^ dr 
and {Dz)'^y is a path integration with boundary conditions imbedded (this is 
marked by the subscript (xy)) and with all possible windings in the Euclidean 
temporal direction. 

One can write it explicitly as 

Here ({k) = z{k) - z{k -l),Ne = t. 



K 
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One can check that in the free case, -B^ = 0, Eq.(12) reduces to well- 
known form of the free propagator 

POO rl^T^ 

= y / e-^''*-^P(^-^)-'P^"^dt-^ (15) 

(2t:Y ^ ' 

dC{m) 



(27r) 

with 



dam) 



Using the Poisson summation formula 

^ ^ exp{ip4n(3) = ^ 5{p4,(5 - 2'Kk) (16) 

^"^ n=0,±l,±2... fe=0,±l,... 

one finally gets the Jackiw-Templeton's form [1,2] 

Note that as expected the propagator (12), (17) corresponds to a sum of 

ghost paths with all ppossible windings around the toriis. The momentum 
integration in (14) asserts that the sum of all infinitesimal "walks" ({m) 
should be equal to the distance {x — y) modulo n windings in the conpactified 
fourth coordinate. 

For the gluon propagator in the backgrounvd gauge we obtain similarly 
to (12) 

{W)Zl = j^ dt{Dz):^e-''^Fix,y) (18) 

where 

^F{x,y) ^ PFPexp{-2ig / F{z{T))dT) expig / Bf^dZf^ (19) 

J Jy 

and the operators PpP are used to order insertions of F on the trajectory of 
the gluon [11,12]. 
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Now we come back to the first term in (11), < (f >b, which can be 
represented with the help of (12) and (18) as 

< >B= J jC{t)d'^x{Dz)'^^e~^[^tr < ^f{x,x) >b - < tr^x,x) >b] 

(20) 

where the sign tr implies summation over Lorentz and color indices. 

In the Appendix 1 we show that Eq.(20) yields for 3^ = the usual result 
of the free gluon gas: 

F,{B = 0) = -T<p{B = 0) = -{N', - 1)^3^ (21) 

As a next example we consider the case of nonzero component -B4, while 
Bi = 0,i = 1,2,3. We also neglect for simplicity the difference between 
and $, i.e. we put F = 0. Actually the interactions of F represent the 
interaction of the gluon spin with the background field [11, 12] and gives rise 
to short-range correlations (in the confining phase it also yields spin-orbit 
interaction, the nonperturbative part of which, the Thomas precession, is 
long range, ~ 1/r) with the notation 

/•/3 . 

< Q >B=< expig / BJz)dz4 >b (22) 
Jo 

one obtains from (20) 

<^>B=TrJ^am(3^^ y: e-^''-'^^-^<cr>B- 

= -TrM r de{l + 2 Y: e-®" < f^" >b) = (23) 

= Trc{-2Vs / ^ < ln{l - e-^^) >b= 

^Tr^<— . >B= Ttc^Vs 2^ — 24 

In the case Q = 1 we come back to the free case, Eq.(21), since 

00 1 4 

E 4 = - (25) 
t,n' 90 ^ > 
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The expression (24) for a constant field B4 yields the well-known result 

Fo{B4 = const) = 5— Trc > j (26) 

^ n=M,... ^ 

where Tvc means trace over color indices. 

One can notice in (26), (24) that in general the presence of field S4(or 
Q 7^ 1) descreases the modulus of Fq, i.e. is not advantageous from the point 
of view of the minimum of Fq. 

Thus < ip >B= —Fo/T in the dcconfined regime represents a loop contri- 
bution of the gluon gas in the background field. Now we turn to the confining 
regime and again calculate < (f >b- It is convenient to choose another point 
y on the loop in addition to the initial point x, as shown in Fig. 3 and to 
write < (/? >B as 

< ^ >B= f mj jl ^ciV^a;(D.)- (D.')- e«--^ (27) 
where we have introduced the area law for the closed contour average 

< ^{x,x) >= exp{-aaS) (28) 

It is important at this point to stress that in the confining regime the total 
contour C'{x,y,x) shown in Fig.4, made by trajectories z{r) and z'{t) in 
(27), should be a closed contour, or in other words, in the representation 

Dz{t) = H dC{m) ——^exp{ip(^C{m) - {x - y) - nl35^^)) 

m=l n l^'^J 

N' j4 / 

Dz'ir') = n ^CK) E 7Ae^P(^P'(EC'K) - {x - y) - n'(55^)) (29) 

m'=l n' K'^'^) 

one should choose 

n^n' (30) 

In this case both trajectories z{t) and z'{t') are winding in the same way on 
the torus and the minimal surface S lies on the surface of the torus, winding 
around it n times. In the case n ^ n' there appear a piece of length \n — n'\ 
where a gluon is winding and contributing an infinite amount to the free 
energy (since we assume that the free energy of an isolated quark or gluon 
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is ininite). Thus for n = n' one can integrate in (27) , (29) over d'^i^x — y) as 
follows (note that X4 — 1/4 changes in the limits [0/9]) 

d'^{x-y)J2exp[-i{p+p'){x-y)-i{p4+p'4){x4-y4 + n(3)] = (27r)V(p+p') 

n 

(31) 

We obtain ior < (p >b 

<^>B=J at)j^V,pG{h,t-h) (32) 
where we have introduced 

G{h,t- h) = ^.e^'^^^UdCrdCRC-^^ (33) 

and Cr = dm) - C'(m), Cij = ^(C(m) + C{m)) 

One can notice that temperature enters only as a factor f3 in (32) and it is 
cancelled in the free energy 

Fo = -T<^>B=-V, J C{t)j^G{h,t-h) (34) 

One can also rewrite G through the two-gluon Green's function 

G{xx, 00, t,t — ti), where gluons start at the point and meet at the point 

x: 

G{ti,t-ti)= J d'^xG{{xx;00;ti,t-ti). (35) 

In the free case it is equal to tI;;^^"^^*. Note however that this result co- 
incides with the loop diagram for T = 0, which is entirely absorbed by the 
subtraction constant in calculating (21). 

In case of confinement one obtains instead 

G = E^n(0)e-^"* (36) 

n 

where fni^), Mn are wavef unction at origin and glueball mass respectively. 

In any case our result for the gluon loop with confinement, Eq.(34), does 
not depend on temperature and should be included in the value of the free 
energy for zero temperature, which as we shall argue below is given by the 
scale anomaly [14] ,i.e. by the gluonic condensate. 
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Note that physically this correspondence is quite understandable, since 
the gluon loop contributimon (34) is expressed through the two-gluon bound 
state at the total momentum zero, i.e. through the condensate of bound 
gluon-gluon pairs (glueballs), and the latter is most generally given by the 
gluonic condensate < > via the scale anomaly expression. 



4 The lowest order quark contribution 

Integrating over quark fields in (3) is done trivially and leads to the following 
additional factor in (7) 

det{m + D{B + a)) = ]-det{m'^ - D'^{B + a)) (37) 

where we have used the symmetry property of eigenvalues of D In the lowest 
approximation we omit in (37) and write the free energy contribution 
Fo(-B) similary to the gluon contribution (9-10) as 



^-Fl{B) = -Undetim' - D\B)) = -^Sp H C(t)^e-*-^+*^^(^) (38) 
1 Z 2 Jo t 

where the sign Sp has the same meaning as in (10) and 



= (L>^7^)2 = DliB) - gF^^a^, = - gEF; a^, = +^(7^7. - 7.7m); 

(39) 

Our aim now is to exploit the FSR to represent (38) in a form of the path in- 
tegral, as it was done for gluons in (14). The equivalent form for quarks must 
implement the antisymmetric boundary conditions pertinent to fermions. It 
is easy to understand that the correct form for quark is 

- -\tr [ at)jd'x-iDz)le-^-'-'Wj2{Cn) (40) 
where Wj^iCn) = PpPAexpig Jc^A^dZf,expg{Y. F) 

(^;=n^ E (-l)"^e-CL.CM-(^-.)-n/^^.. (41) 

=1 (.^Tit; n=0,±l,±2,... l^^J 



m= 
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One can easily check that in the case 3^ = one is recovering the well- 
known expression for the free quark gas 

F^ifree quark) = -^-^N^V^T^ ■ Uf, (42) 

where Uf is the number of flavors. The derivation of (42) starting from 
the path-integral form (40) is done similarly to the gluon case given in the 
Appendix 1. 

The loop Cn in (40) corresponds to n windings in the fourth direction. 
Above the deconfinement transition temperature Tc one can visualize in (40) 
the appearance of the factor 

= exp ig / B4(z)dz4 (43) 
Jo 

(Note the absence of the dash sign in (43) as compared to (22) implying that 
in (43) B4 is taken in the fundamental representation) . For the constant field 
B4 and Bi — 0,i — 1,2,3 one obtains 

<^>=->.E^^^(-ir' (44) 

^ n=l ''' 

This result coincides with the obtained in the literature [13]. 



5 Higher order contributions 

We unify under this title two different types of contributions: i) higher order 
cumulants from the cluster expansion (11) for the determinantal term (9); ii) 
terms of the higher order in g, appearing in the perturbative expansion (7). 

We discuss first the terms i). A generic case is given by the quadratic 
cumulant <^ ip^ which can be written as (we neglect the gluon spin 
terms) 

mT at')^d'xd'x'e-^-^\Dz)-^{Dz')-,^,x (45) 
Jo t Jo V 

X < tr^{x,x)tr^{x' ,x') > 
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where $ is defined in (13). The cumulant in the integral (45) corresponds to 
the interaction between two gluon loops, each loop being a gauge-invariant 
object. Hence the interaction of the loops, given by the cumulant, is sup- 
pressed as l/Nl both in the confining and nonconfining regimes(as compared 
to the contribution < ip >b in (11)). 

In the confining regime two gluons form a bound state - a glucball. One 
can choose the initial glueball state connecting the loops {xx) and {x'x') by 
the stright line at some fixed x^ — x'^, and then (45) can be rewritten as 

« »= / at)j J at')^Vip\-^-^'-^^'iDR):^iDr)-, (46) 

where R and r are central and relative coordinates of two gluons defined in 
accordance with (33). The resulting integrals in (46) are similar to those in 
(32) with one essential difference; in (32) the integration over all d{x — y) as 
in (31) is performed leading to the total momentum of two gluons = 0. 
In contrast to that in (46) the initial and final central point R is the same 
and therefore all values of total momentum are possible as in the one-gluon 
case, see (15) and Appendix 1. 

As a result one obtains the temperature dependence and the final result 
(see Appendix B for details) is the free energy of the glueball state [15] 

<W»-Y.'-^^^-- (47) 

k 

where is the glucball mass in the glucball state k. Note that (47) is 0(iV°) 
since it is a color singlet contribution, in contrast to the gluon contribution 
(21). 

In the same way one can consider the next irreducible average, <^ 3>. 
The contribution is again a color singlet, corresponds to the three-gluon 
glueball in the confining regime and is of the order 0(A''°). 

The same features are common to higher order cumulants. 

In a similar manner one can study the quark contributions to ^ (^^ 
In this case one has to replace in (45) the adjoint parallel transporters $ by 
the fundamental ones and use the fermion boundary conditions in {Dz)'^y as 
in (41). In the confining regime one obtains the meson contribution to the 
free energy which has the same form as in (47). It is of the order of (A^^)- 
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The largest contribution at small temperatures, T <^Tc comes from the pion 
gas, where neglecting the pion mass we have 



Fo{pion gas) = -— FgT^ 



(48) 



One should note, that in the set i) introduced above the gluons interact 
whith each other via nonpcrturbative field 5^, which can be exemplified by 
the correlators <ti F^i,{B)F\f^{B) ^ etc. In contrast to that in the set ii) all 
perturbative exchanges are taken into account. These are proportional to g'^ 
and in the usual treatment g — g{T) is small at large temperature while it 
may diverge as T ^ A [1] 



In the strong background field however one has an additional scale pa- 
rameter < trF^^{B) >= B which defines dynamics both at small and large 
distances and temperatures. One can show that when B ^ \^ this scale 
defines the value of the coupling constant g{B,T) and g{B,T) remains small 
even when T tends to zero [7]. Therefore one may hope that the perturbative 
terms (set ii)) which yield relatively small contribution (~ 10%) at T ~ 3Tc, 
remain of the same order at all temperatures. 



6 Calculation of the deconfiniment tempera- 
ture in the leading order of the l/N^ expan- 



In this section we estimate the deconfinement temperature keeping only lead- 
ing terms in both phases - confining and deconfining. Results given below 
have been published in a short form in [8]. We start with the zero tem- 
perature and note that the NP energy density e is connected via the scale 
anomaly [14] to the gluonic condensate [16] (we neglect the quark masses 
since only light quarks are important at low temperatures) 



AT) 



247r2 



(49) 



{llNc-2Nf)lnT/A' 



sion 



e — 



P{as) 



(50) 



16a 



c 



327r 



13 



One can associate with e the zero-temperature Umit of the free (Helmholtz) 
energy F(T = 0) and write F at nonzero temperature as 



F^eV^ + f{T) 



(51) 



where f{T) is to be computed using the formulas (8) and subsequent ones. 
Since f{T) is regularized by subtracting the zero-temperature limit /(T = 0) 
and all divergencies are present in the latter term, Eq.(51) is a definition of 
the regularized limit F{T = 0) where only NP contribution in the gluonic 
condensate is present. A check of correctness of our normahzation oi F{T — 
0) is given by the computation of the leading contribution in Section 3. 

There the perturbative contribution (21) goes to zero as T" — > in accor- 
dance with our prescription that perturbative contributions are normalized 
to vanish at T — > 0. If however one takes into account the confining back- 
ground, one gets for T ^ a nonzero contribution (34), which is exactly due 
to a pair of gluons connected by the adjoint string (28) and with zero total 
momentum (see (31)). 

This is exactly the contribution which one expects from a gluon pair 
forming a condensate. 

The low-temperature phase consists of the gluonic condensate (50), glue- 
balls (47), mesons (we keep only the pion (48)) and their interacting con- 
glomerates: 



Here the first term is the leading in N^, it grows as [17], while the ideal 
gas of glueballs and mesons contributes 0(iV°). Note that the interaction 
between white objects is suppressed at N(. ^ oo [18] and is presented in (52) 
by the last term. 

The high-temperature phase is belived to be the gas of quarks and gluons, 
interacting perturbatively [1,2]. We shall now argue that i) there is at T > 
another component of the phase, namely a part of the gluonic condensate, 
which has not evaporated during the phase transition; ii) there is a NP 
interaction between quarks and gluons which is important in some situations. 

To start we remind that the NP interaction is governed by the correlators 



Fi^^^eVs-TY^ 



K 



87r3/2 



- ^^3T^ + 0(1/Ar,) (52) 



[6] 



G^^,xa =< tr[F^^{x)^{x,y)Fxa{y)^{y, x)] > 



(53) 
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which at T = contain two independent Lorentz structures with scalar 
coefficients D{x — y) and Di{x — y) [5,6] 

1 d 

Gpi^,au{u) = {SpaSi^u-SpuSHa)D{u) + -[—{{u)aSij,u-{u)Ji^a) + {P'^ ^ IIiy)]Di{u) 

(54) 

At T" > the Lorentz invariance is broken and one obtains two independent 
correlators- for electric and magnetic fields respectively [19] 

Glf = G,4,fe4 = S^.D^'^Xu) + Ii..]d[''\u) (55) 

Gjf^ = Gem,pq^eiem^ekpq = 6ikD^ (u) + ^[...]Df{u) (56) 

For the gluon condensate one must put ti = in (54-55) with the result 

< irF^,(0)F^,(0) >= D^{0) + Df{0) + L'^(O) + Df{0) (57) 

Out of four functions in (56) only one, D^, is connected to the confinement, 
since the string tension a is expressed through it [5,6] 

(7 = ^ y D^{^ul + ul)duidu4 + ... (58) 

where dots imply contributions of higher-order cumulants. 

Note that the integral (58) is confined to the region ^Ju\-\-ul -< Tg, since 
falls off exponentially outside. The value of Tg - the gluonic correlation 
length - is the order of 0, 2/m [20]. 

Therefore for the temperatures T ^ T^. — 200MeV the periodicity in 
the fourth direction does not change a significantly, and the deconfinement 
occurs when disappears at some T ~ Tc together with a. There is no 
reason why Df, D^, Df disappear at T = Tc or at T > T^. Moreover, there 
is evidence [21] that D^, which defines the area law of the spacial Wilson 
loops, should be nonzero at T > T^. Additional evidence comes from the 
hadronic screening lengths [22], also implying the confining dynamics in the 
spacial directions (we refer the reader to [9,19] for more discussion of this 
point). 

Thus we can assume that only disappears at T = Tg. How gluonic 
condensate changes with temperature? The answer comes from lattice calcu- 
lations [23] and shows that < trF^j^(O) > is rather stable and only a part of 
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it disappears near T = Tc- This is understandable from the theoretical point 
of view, since the natural scale of < trF^(O) > is given by the lowest glueball 
mass approximately equal to the mass of the dilaton [24] which is of the order 
of IGev, and therefore < irF^(O) > should not strongly change in the range 

^ T ^ Tc . Thus we come to the following model for high-temperature 
phase: 

FugH = (1 - n)eV, - {Nl - - ^^cV,T%f (59) 

where 

Since Di enters the NP tensor force of quarkonia [11] one can make some 
estimate oi Di/ D : Di/D ^ 1; the lattice calculations [20] at T = yield 
%^ ^ In > 0.1 fm. 

Therefore it seems reasonable to neglect Df and Df ( note that Df — 
Df ^ Di ai T ^ 0) and having in mind that ^ Sit T ^ 0, one has 
an estimate 

V = 1/2 (61) 

As another option with D^{0) = -Df (0) = i:>^(0) = -Df (0), one can choose 
7] = 1/4. 

We have neglected in (59) i) higher order terms in g"-, ii) higher order 
cumulants <S ^ with n >z 2 and iii) the influence of the term Q (24) 
and (43). As we discussed above, we believe that i) ag is not large, ^ 0.5 
at all temperatures ii) higher cumulants are suppessed as 0(1/A^'^), iii) from 
the structure of (24) and (43) one can notice that Q ^ 1 tends to decrease 

1 F I and since in the equilibrium the phase chooses the minimal F (maximal 
I F I) it will choose by that Q — 1. Note also, that studies of F for constant 
-B4 [25] seem to support our conclusion (in lowest orders graphs). 

Of course our argument here is very qualitative and can be used only for 
a rough estimate of T^. A more detailed study of the role of Q is now in 
progress [26]. We now can plot both curves Fhigh (59), Fiow (52) or rather 
the pressure P = — F/V3 as a function of temperature and use again the 
principle of minimal F (maximal P). One can see in Fig.4 that there exists 
indeed some critical temperature T = Tc and the system prefers to be in the 
low phase at T < Tc and in the high phase at T > Tc. The transition is of 
the first order, and Tc can be computed equalizing (59) and (52) 
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In the leading order 0{N^) one obtains: 



We note that Tc is 0(A^°) which is reasonable since strong interactions (ex- 
cept for baryons) have finite limit both in masses and the range of interaction 
when Nc ^ oo. With the standard value [16]. 

G2 = — <G^>= 0M2GeV* (63) 

TT 

we have for rj — 1/2 

Te = (0.196r;G2)^/^ = O.lSbGev (64) 

However number of flavours Uf — gluonic gluonic condensate should be 2-3 
times larger [16], which yields Tc = 0.22 — 0.24:GeV and agrees well with the 
lattice calculations for gluodynamics [9]. 

In the next order O(A'c) we take into account the quark contribution and 
still disregard all meson and glueball contributions (as well as nonideal gas 
corrections for gluons and quarks). We obtain 

11 j\T as<G'^> 

For the choice (61) and Uf — 2 we have T^. = O.lbGeV while for Uf — A 
we have = 0.134G'ef , which agrees well with recent lattice calculations 
[9,10] T^^^-^ (lattice) =0.14aGey and 0.131^6^ respectively 



7 Conclusion 

We have started in the paper the background fleld formalism for T > with 
the purpose to take into account the NP contributions in a systematic way. 

The perturbation theory was formulated where expansion is of two dif- 
ferent kinds: the cluster expansion for nonperturbativc interaction and the 
usual expansion in powers of the coupling constant for perturbative quarks 
and gluons in the strong NP background. 
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In this first paper we concentrated on the lowest order terms: i) the 
gluonic condensate in the scale anomaly term ii) gluonic and quark loop 
diagrams in the background field. 

We have shown in the previous Section that already a rough approxim- 
tion, where only terms i) and ii) are kept, yields a reasonable picture of the 
deconfining phase transition where a) the order of the transition b) the nu- 
merical value of Tc c) existence of non-perturbative effects at T > Tc, are 
correctly predicted. We also demonstrated that the terms i) and ii) are lead- 
ing in the large Nc hmit and the resulting value of Tc doesn't depend on Nc 
in this limit. 

We plan to expand the formalism to include higher order terms and check 
the stability of our predictions. First of all, the role of gluonic (adjoint) 
Polyakov fine Q is being clarified [26]. This role is especially important for 
the question of the latent heat, which comes out too large with Q = 1 as 
compared to lattice calculations. 

Another important possibility of the presented formalism is the study of 
infrared singularities of the free energy for T > Tc connected to difficulty 
with the magnetic mass [3,13]. 

In the vaccum background field of QCD the confinement BiT > Tc is still 
present in the special planes [19] which cures the infrared divergencies and 
can solve in principle the problem of magnetic mass. This will be discussed 
in detail in a subsequent paper. Note added: after the original version of this 
manuscript has been published as a preprint [27], the author became aware of 
a similar estimate of Tc in [28] , where it was assumed that no nonperturbative 
configurations are present for T > Tc in the high temperature phase. This 
corresponds to r] = 1 and results in ~ 25% higher values of Tc as compared 
to ours. 

However in the pictures of [28] it is difficult to explain lattice results of 
[21-22] which require nonperturbative magnetic configuration for T > Tc. 

The author is grateful for discussions to the members of the theoretical 
seminar of ITEP. A part of this work has been done while the author was 
a guest of the Max-Planck Institut fiir Kernphysik in Heidelberg. It is a 
pleasure for him to thank Professor H.A.Weidenmueller and all the staff 
of the Insitute for a kind hospitahty. Fruitful discussions with H.G.Dosch, 
H.Leutwyler, P.Minkowski and H.D.Pirner are gratefully acknowleged. 
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APPENDIX 



Calculation of the free energy of a noninteracting gluon loop via 
the path integral, Eq.(20). 

The square brackets in (20) yield ■ A - 1){N^^ - 1) = {N^ - 1). The 
integral {Dz)"^^ can be calculated as the integral in dC,{m) (see Eq.(15)) and 
using also (16) we get 



Vsd^p^ fl3^p' de 



2 

(A1.2) 



k=o,±i,.A^^r ^1 e2 + (27rfc)2 
Using a relation [1] 



we obtain 

(f f Vsd^p rm , 2 



- y + 2/n(l - e-'^^) + const] = (Al.5) 
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FIGURE CAPTIONS 



Fig.l. The gluon loop diagram in the background field (the first term 
onther.h.s. of Eq.(lO)). 

Fig. 2. The ghost loop diagram in the background field (the second 
term on the r.h.s. of Eq. (10)). 

Fig. 3. The gluon loop with the confining film as a result of averaging 
over background field - Eq.(27). 

Fig. 4. Paths of integration z{t) and z'{t) in Eq.(27) in the in the 
confining regime. 

Fig. 5. The pressure Piow (corresponds to Fiow Eq.(52)) as a function 
of temperature - dash- dotted line. The pressure Phigh corresponding to 
Fhigh, eq.(59)) with the quark contribution -solid line, without it - dashed 
line. Transition temperatures Tc{q) and Tc obtained with and without quark 
contribution respectively. 
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